Contact lines for fluid surface adhesion 
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When a fluid surface adheres to a substrate, the location of the contact line adjusts in order to 
minimize the overall energy. This adhesion balance implies boundary conditions which depend on 
the characteristic surface deformation energies. We develop a general geometrical framework within 
which these conditions can be systematically derived. We treat both adhesion to a rigid substrate 
as well as adhesion between two fluid surfaces, and illustrate our general results for several impor- 
tant Hamiltonians involving both curvature and curvature gradients. Some of these have previously 
been studied using very different techniques, others are to our knowledge new. What becomes clear 
in our approach is that, except for capillary phenomena, these boundary conditions are not the 
manifestation of a local force balance, even if the concept of surface stress is properly generalized. 
Hamiltonians containing higher order surface derivatives are not just sensitive to boundary transla- 
tions but also notice changes in slope or even curvature. Both the necessity and the functional form 
of the corresponding additional contributions follow readily from our treatment. 

PACS numbers: 87.16.Dg, 68.03.Cd, 02.40.Hw 
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I. INTRODUCTION 

There exist a number of physical systems whose ener- 
getics is fully described by a surface Hamiltonian. The 
easiest and best known examples involve capillary phe- 
nomena p], Q , where the shape of a liquid- fluid interface 
(such as a sessile water droplet or the shape of a fluid 
meniscus) is determined by minimizing its area. The 
same physics governs the behavior of soap films. Higher 
order surface properties, notably its curvature, play a 
role in the description of fluid lipid membranes or mi- 
croemulsions [3|, |J] , and even higher derivatives have been 
implicated in the occurrence of certain corrugated mem- 
brane phases Q . In all these cases the shape of the sur- 
face follows from minimizing the surface Hamiltonian, a 
variational problem. The corresponding Euler-Lagrange 
differential equations are known as the shape equations. 

However, such surfaces are generally not isolated but 
rather in contact with something else. Water droplets 
or lipid membrane vesicles may rest on a substrate, and 
this generally influences their shape quite strongly. For 
instance, water droplets on hydrophilic substrates (e.g. 
clean glass) resemble flat contact lenses, while on very hy- 
drophobic substrates (e.g. teflon) they are almost com- 
pletely spherical. When gravity can be neglected @ , the 
shape equation dictates a constant mean curvature sur- 
face in both cases (in fact, a spherical cap), but the con- 
tact angle at the three-phase line where water and sub- 
strate meet is different for the two different substrates. 

In the majority of cases the spatial extension of the 
surface being studied exceeds the range of interaction 
between it and some substrate by a large amount. For 
instance, van der Waals forces, hydrophobic interactions, 
or (screened) electrostatic forces typically extend over 



several nanometers, while the extensions of vesicles or 
droplets can be microns or even millimeters. Under these 
conditions the interaction is well approximated by a con- 
tact energy, i. e., an energy per unit area, w, liberated 
when the surface makes contact with the substrate. It is 
this adhesion energy, together with the energy parame- 
ters characterizing the contacting surfaces, which deter- 
mines the boundary conditions holding at the contact 
line. In the case of capillary phenomena for example the 
ratio between adhesion energy w and surface tension a 
determines the contact angle i9 between liquid and sub- 
strate surface by means of the well-known Young-Dupre 
equation [l|, [|J 
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The "standard" derivation of this result involves a bal- 
ance of tangential forces at the contact line. Yet, despite 
being very intuitive, this requirement of surface stress 
balance does not yield the correct condition for more 
complicated surface Hamiltonians, even if the concept 
of surface stress is generalized properly. Higher order 
Hamiltonians give rise to additional energy contributions 
when the contact line is varied. It is the purpose of the 
present article to show how these contributions can be ac- 
counted for in a systematic and parametrization-free way, 
and without assuming any additional symmetries (such 
as axisymmetry or translation symmetry along the con- 
tact line). We study adhesion to rigid substrates as well 
as to dcformable surfaces also characterized by a surface 
Hamiltonian. Our presentation generalizes and strongly 
simplifies the analysis previously given in Ref. Q- 
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FIG. 1: Illustration of the geometry of surface adhesion. Per- 
pendicular to the contact line C two vectors exist, I and I, 
which are tangential to surface S and substrate 5, respec- 
tively. 



II. MATHEMATICAL SETUP 

A. Differential geometry 

In order to describe the adhering surfaces in a 
parametrization-free way, we use a covariant differential 
geometric language. Our notation is essentially standard 
and will follow the one used in Refs. 0, H, H E3, EH- 
Briefly, a surface S is described by the embedding func- 
tion X(£ x ,£ 2 ) S M 3 , where the £ a (a S {1,2}) are a 
suitable set of local coordinates on the surface. This in- 
duces two local tangent vectors e a = dX/d£ a = d a X 
and a normal vector satisfying n ■ e a = and n 2 = 1. 
Furthermore, the two fundamental forms of the surface 
are needed, namely (i) the metric tensor g ab = e a • ej, 
and (ii) the extrinsic curvature tensor K ab = e a ■ V&n = 
— n • V a eb. The symbol V a is the metric-compatible co- 
variant derivative. The trace of the extrinsic curvature 
tensor will be denoted by K = = K ab g ab , which for a 
sphere of radius r with outward pointing normal vector 
is positive and has the value K = 2/r. As usual, indices 
are lowered or raised with the metric or its inverse, re- 
spectively, and a repeated index i (one up, one down) 
implies a summation over i = 1,2. More b ackg round on 
differential geometry can be found in Refs. jl2l T13l. [LH]. 



B. Geometry at the contact line 

By "contact line" we will denote the curve C at which 
the surface detaches from the substrate. Its local direc- 
tion is given by the tangent vector t = t a e a (see Fig. [TJ, 
which is tangential to C, the surface 5, and the substrate 
surface S_ (which itself might also be deformable). Here 
and in what follows (with the exception of Sec. IIVB|) we 
will use underlining in order to indicate quantities refer- 
ring to the substrate. 

Perpendicular to C we can define local normal vectors 
which are either tangential to S or <S, namely I — l a e a 
and I — / a e a , respectively (see Fig.[T|). Also, we will have 



two surface normals n and n. If the surface contacts the 
substrate at zero contact angle, we will have e a = e a , 
n = n, and I — I there; however, their derivatives per- 
pendicular to the contact line need not coincide, since the 
curvatures of surface and substrate generally need not be 
identical. In fact, the values of perpendicular and par- 
allel components of these curvatures (and possibly their 
higher derivatives) will be among the primary focus of 
this paper. They will be denoted by 





= l a l b K ab = 


—n 


v ± z , 




= t a t b K ab = 


—n 


V||t , 




= l a t b K ab = 


—n 


V±f = -n 



(2c) 

where we also introduced the two directional surface 
derivatives perpendicular and parallel to C, 

V ± =rV a and V||=t Q V a . (3) 

Notice that we can analogously define K ± , Ku, and K_±\\ ■ 

C. Hamiltonian 

In the present work we will exclusively study surfaces 
whose energy H is given by a surface integral over a scalar 
energy density 7i, which is constructed from the local 
surface geometry: 

H = [ dA H(g ab , K ab , V a K bc , . . .) . (4) 
Js 

The integral extends over the entire surface 5, and the 
area element d^4 is given by y^d^d^ 2 , where g = 
det(g ab ) is the metric determinant. For those parts of 
the surface which adhere to a substrate we will assume 
an additional adhesion energy density 

Wadhe S ion = -W^ 1 ,?) <0 , (5) 

which may in general be a function of position. 

III. DETERMINING THE BOUNDARY 
CONDITIONS 

A. Continuity considerations 

As we will see, the adhesion balance between surface 
and substrate will result in a discontinuous change of 
some surface property across the contact line. However, 
the form of the energy density restricts which quantities 
can be discontinuous, since it needs to remain integrable. 

Most obviously, the shape itself has to be continuous. 
Yet, already its first derivative may display a jump, as it 
does in the case of capillary adhesion. The energy density 
is given simply by 

^capillary = CT , (6) 
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and a kink in the surface at the contact line, i. e. a finite 
contact angle, is not associated with an extra energy. 

For curvature elastic surfaces the situation is different. 
There, the energy density of the surface is given by the 
well known expression 0, 0] 



Hbcnd = - Kof 



(7) 



where n is the bending modulus (with units of energy) 
and Ko describes a spontaneous curvature of the elastic 
surface [l5j]. A kink in the surface at the contact line 
implies a ^-singularity in the curvature, whose square is 
non-integrable. Hence, the surface needs to be differen- 
tiable across the contact line, and the distinction between 
surface- and substrate tangents and normals drawn in 
Sec. Ill Bl becomes unnecessary. Moreover, a quick glance 
at Eqns. (|2bl2cp shows that both Kn and the off-diagonal 
curvature K±u is expressible as a tangential derivative 
along the curve C of a quantity continuous across C, hence 
both these curvatures will also be continuous. It is only 
the perpendicular curvature component K± which might 
possess a discontinuity, and indeed we will see that it 
does. 

Finally, even higher order derivatives might occur in 
surface Hamiltonians. For instance, Goetz and Hclfrich 
have studied a curvature-gradient term of the form 



B. Contact line variation 

For an adhering surface the total energy is station- 
ary with respect to variations of the contact line along 
the substrate. Such a variation contributes twofold to 
the Hamiltonian: Assume that locally the contact line is 
moved such that a bit of surface unbinds from the sub- 
strate. This removes its corresponding binding energy, as 
well as any elastic energy associated with the constraint 
of conforming to the substrate, and thus gives rise to 
an energy change <5-fft>ound- On the other hand, the un- 
bound part of the surface acquires at the contact line a 
new boundary strip which implies also a change 5Hf lee in 
its elastic energy. The boundary condition at the contact 
line then follows from the stationarity condition 



SH cl = SH h 



ound 



SH flcc = 



(10) 



In the case of adhesion to a rigid substrate the bound con- 
tribution involves the variation along a surface of known 
shape. The corresponding term is thus conceptually very 
different from a deformable substrate or even the free 
variation, because in both these cases the local shape of 
the surface is not known. Below we will see how these dif- 
ferences manifest themselves when computing the bound- 
ary terms. 



n. 



grad = ^K V (V a K)(V a K) = ^K V (VKf , (8) 



which in a generalized higher-curvature Hamiltonian pre- 
vents the occurrence of infinitely sharp curvature changes 
Q. In this case it is obvious that all curvature compo- 
nents have to be continuous along the substrate, since 
otherwise again a squared <5-singularity results. More- 
over, most of the first order directional derivatives are au- 
tomatically continuous: The parallel ones, V||, again dif- 
ferentiate quantities along C which are continuous across 
C and thus are itself continuous. For the perpendicular 
ones it turns out that Vx-K) and V_l-Kj_|i are continu- 
ous, while W±K± is not. This is intuitively reasonable, 
since every term involving a "||" features at least one 
less derivative across the contact line and thus cannot 
jump. A rigorous proof is however a bit more involved. 
One may for instance proceed like this: Start with 
the contracted Codazzi-Mainardi equation [13, [IH G3 
V a K ab — V b K = and project onto lb- By decomposing 
the resulting identity into the local (I, t) frame, it can be 
cast in the form 



V 



X,| =V||Ai||- 



■{K±-K ll )t-'V\\l- 



-2if_L||Z-Vj 



t. (9) 



Since every term on the right hand side is continuous 
across C (recall that derivatives of tangent vectors are es- 
sentially curvatures), \7±K\\ must be continuous as well. 
By projecting the contracted Codazzi-Mainardi equation 
on tj, instead, one can show continuity for V ' _\_K±_\\. 



1. The bound variation 

For definiteness, let the normal I to the contact line C 
be directed towards the adhering portion of the surface 
(see Fig. [T]). A local infinitesimal normal displacement 
e of the contact line along a rigid substrate thus implies 
the following obvious change in the bound part of the 
surface: 



5H h 



ds (n 



w)e(s) 



(11) 



The underlining of TL should again indicate that it is eval- 
uated with geometric surface scalars (such as for instance 
curvatures) pertaining to the substrate. If the substrate 
is flexible, the w term remains, but the change in elastic 
energy will instead be taken care of by an additional free 
boundary variation. 



2. The free variation 

The change in energy due to the addition or removal of 
unbound parts to the boundaries of the surface is identical 
to the boundary terms in the variation of the free surface. 
In Ref. 9] it has been shown that for Hamiltonians up to 
curvature order these terms are given by [l6j 



SHf ICC — 



- J ds l a 



f a -6X 



H ab n 



Se b 



(12) 
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Here, f a is the surface stress tensor, given by 

f a = {T ab - H ac K b c )e b - {V b H ab )n , (13) 
and we have also defined 

sn 



H 



ab 



T 



ab 



and 

dKab 

V9 $9ab 



(14a) 
(14b) 



Finally, 5X and 5e b denote the change of contact line 
position and the associated change in the slope of the 
tangent vectors, respectively. Notice that the latter term 
is only relevant if TC ab ^ 0. 



IV. SPECIFIC EXAMPLES 

In this section we will illustrate the above formalism 
by applying it to several important situations and surface 
Hamiltonians. In Sec. IIV Al we first treat the problem of 
adhesion to a rigid substrate. We will see how known 
results (the Young-Dupre equation and the contact cur- 
vature condition for Helfrich-membranes) follow with re- 
markable ease and can be extended just as quickly to 
new Hamiltonians. In Sec. IIV 51 we look at the boundary 
conditions involving adhesion to deformable substrates. 
Specifically, in IIV B ll we look at the triple line between 
three tension surfaces, and in IIV B 2l we study the adhe- 
sion of two vesicles. 

A central ingredient in all this will be the knowledge of 
the two tensors H ab and T ab defined in Eqns. (|14all4bj) . 
While their determination is not particularly involved, 
these calculations have been performed previously by us 
d, [H, Ell EH and we will thus simply reuse the results 
here. 



A. Adhesion to a rigid substrate 

Since the variation of the contact line has to proceed 
along the substrate, we must have 



SX = el 



(15) 



No component in t direction is necessary, since for fluid 
surfaces this would merely amount to a reparametrization 
of C. Notice that (fl~5|) is nothing but the Lie derivative of 
X along the substrate, since C e iX = el_ a W a X = e{ a e Q = 
el. This property holds generally, and we will make use 
of it later. 

The normal component of the change in the surface 
tangent vectors e b only contributes if Ti ab ^ 0, i. e., if 
curvature terms enter the Hamiltonian. We will assume 
that they do it in such a way that differentiability of the 
surfaces is implied (see Sec. lIII A]) , so that no distinction 
needs to be drawn between normal and tangent vectors 
of substrate and adhering surface. We then find 



where in the last step the equation of Weingarten [12|, 
EJ, EH has been used; this is again the Lie derivative 
along the substrate • Notice that there still remains a 
distinction between curvatures of substrate and surface; 
hence, the derivative of the tangent vectors resulting from 
a variation along the substrate yields the substrate cur- 
vature and not the free surface curvature. 



1. Capillary surfaces 

In this case the energy density is given by Eqn. 
and as we have seen in Sec. IIII Al we will expect a dis- 
continuity in slope at C. The bound variation is 

Abound = - J ds (cr - w)e(s) . (17) 

For this Hamiltonian we have [1, ©, ES EH 'H ab = and 
/" = —ageb = — ere a , and therefore 

l a f a -5X = -al a e a -el = -osll_ = -ae cos(tt-i9) , (18) 

where "& is the angle between capillary surface and sub- 
strate - in other words, the contact angle (see Fig. [l}. 
Equation (fT0|) thus specializes to 



SH cl = - / ds [cr - w + a cos #] e(s) . (19) 
Jc 

Since e(s) is arbitrary, the term in square brackets must 
vanish - which gives the Young-Dupre equation JT]). 



2. Helfrich Hamiltonian 

Let us now look at the energy density ([7]) which de- 
scribes the continuum behavior of (tensionless) fluid lipid 
bilayers. In this case 7i ab = n(K — Ko)g ab and for the 
stress tensor we have B S Ed, m 



U a ■ l = ~-k(K - K j 2 + k{K ~ K )K ± . (20) 

Together with the expressions for the contact line po- 
sition and tangent variation from Eqns. (|15ll6p we find 



SH cl = -J ds ^k(K-Ko) 2 -w 



-k{K - K ) 2 + k(K - K )K A 



-k(K-K )K x \ s(s 



(21a) 



n ■ Set = n ■ Vft<5X = el c n ■ V&e c 



-el c K 



he 



(16) 



= - J^&s ^l K (K ± -K ± ) 2 -w\ e(s) , (21b) 

where in the second step we made use of the continuity 
condition K\\ = K_u discussed in Sec. IIII Al Equation 
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(121b[) implies a discontinuity in the perpendicular curva- 
ture K± as the appropriate adhesion boundary condition: 



(22) 



K±-K, — \ — 



The correct sign after taking the square root follows from 
the fact that the detaching surface must not penetrate 
the substrate; unfortunately this depends on ones specific 
choice of the surface normal vectors. 

Quite remarkably, this boundary condition depends 
neither on the spontaneous curvature Kq, nor on the local 
parallel curvature K\ | . It would also remain unaffected 
if the bilayer were under a finite tension a . Formally, it 
is easily seen to cancel; physically, the reason is that the 
jump we would expect from the Young-Dupre equation 
(fTj) cannot materialize since the curvature terms in the 
energy density enforce differentiability of the profile at C. 

Equation (f2"2l coincides with the result given previ- 
ously in Ref. nfl. Its axisymmetric version was first 
quoted in Ref. 18], and its specialization to a straight 
contact line can be found in Ref. [l9(. We want to stress 
that the H a n-Seb term in Eqn. (fT2"|) . which is responsible 
for the third line in Eqn. (|2f aft , was crucial in obtaining 
equation l|22[). Leaving it out - i.e., only treating the 
problem as a stress balance - will not result in the cor- 
rect boundary condition. The only exception (treated 
via stress-balance in Ref. Q) is the special case of a flat 
substrate, in which case Se^ = and the missing contri- 
bution vanishes anyway. 



3. General Hamiltonians involving K 

The analysis in the previous section readily extends to 
Hamiltonians which are expressible in terms of the total 
extrinsic curvature, H = H(K). One finds H ab = H'g ab , 
where the prime denotes differentiation with respect to 
K, and the tangential stress tensor contribution is 



Lf a ■i = -n + n'K± 



(23) 



The remarkably simple general boundary condition then 
reads [H 



(H-H)-H'{K X -K A 



(24) 



Notice that unlike the form it takes for the quadratic 
Hamiltonian discussed in Sec. II V A~2l this condition gen- 
erally involves the parallel curvature K\\ on the contact 
line. 



4- Curvature gradients 

Finally, we want to study surfaces described by the 
Hamiltonian density ([5]), i. e., involving curvature gradi- 
ents. From earlier work fiol [ill we know that 
n ab 

l a f a ■ I — K V 



Kv (AK)g ab 



and 



(25a) 



(ViT) - Kj_AK .(25b) 



However, simplemindedly inserting these expressions 
into the formulas we have used so far does not give 
the correct result. Here is why: When calculating Tt ab 
and T ab to obtain Eqns. (gSJ), we varied the Hamilto- 
nian ([8]) with respect to K a b and g a b, respectively, and 
identified the bulk terms (see Sec. 4 in the Appendix of 
Ref. (llj|). These variations also leave boundary terms, 
since the curvature K = g ab K a b appears differentiated. 
For the purpose of identifying Ti ab and T ab they are ir- 
relevant, but they evidently matter now that we are in- 
terested in the total energy change upon displacing the 
boundary. Moreover, the auxiliary framework introduced 
in Ref. [9j teaches us that the SK a b and Sg a b variations 
are indeed independent from the SX and Se a terms al- 
ready included in Eqn. (|12|l . so they can be simply added 
there. If we go through the calculation, we see that 
the two tensor variations can be nicely combined into 
a single scalar one, since they occur in the combination 
g ab SK a b + K ab 8g ab = SK, and we end up with the addi- 
tional boundary contribution 



6H, 



frcc.grad 



k v / ds l a {V a K)5K . 
Jc 



(26) 



This is of course exactly the boundary term we would 
expect for the variation of a Hamiltonian density whose 
functional form is the square of the gradient of a scalar, 
so everything is consistent. 

To evaluate the right hand side of Eqn. (|2l)|) , note that 
the variation is once more given by the Lie-derivative 
along the substrate. Since if is a scalar, we obtain the 
simple expression 



SK = C e iK = el a \7 a K = eV±K 



(27) 



Together with Eqns. (|25|) we then obtain the total contact 
line variation as 



SH C 



ds i U v (VKf 



-Ky 



(V±K) 2 - -(WO 2 - K ± AK 



-k v {W±K)[W±K)^ e(s) (28a) 

ds |i/s v (Vxif± - V±K ± ) 2 - wj ep^b,) 

where in the last step we used the continuity of curvatures 
and their V||-derivatives (as discussed in Sec. IIII A|) as 
well as the decomposition (VK) 2 = (V±K) 2 + (V||1T) 2 . 
The boundary condition following from this specifies a 
jump in the perpendicular derivative of the perpendicular 
curvature 



V±(K A 




(29) 



The similarity with Eqn. (|22p is quite striking, and one 
might surmise a pattern that would be followed by even 
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FIG. 2: Illustration of the geometry at the contact line be- 
tween two adhering vesicles. 



higher derivative theories. Notice, however, that the 
terms entering the derivation of Eqn. (|29|) are quite 
different and that the additional term stemming from 
Eqn. (|26[) , which is absent in the simple curvature square 
case, is essential. 



B. Adhesion to deformable surfaces 

Compared to the previous section, there are two key 
differences if the substrate is not rigid. First, the absence 
of a known substrate shape along which a certain amount 
of deformation energy is to be paid removes the term 
involving "H m the bound variation (|11|) . Second, for 
the same reason the contact line variation is no longer 
restricted to proceed along a substrate and will thus be 
of the more general form 



5X 



I 



- s n n 



(30) 



The corresponding tangent vector variation, which occurs 
if TL ab 7^ 0, then leaves a term 



Se b = n ■ V b (e±l + e n n) 



-KbJ c £± + V fc e„ . (31) 



Note that K_ bc plays a different role here than previ- 
ously. It no longer describes the curvature of the evi- 
dently nonexistent substrate. Rather, the tangential vari- 
ation may proceed locally along a fictitious surface which 
is tangential to the other three surfaces that meet at the 
contact line. Encoding higher order derivative informa- 
tion necessary here, K_ bc describes the curvature of that 
fictitious surface, and K ± is the component perpendicu- 
lar to C (see Fig. [5]). This surface is of course not unique, 
and thus K \ is arbitrary - just as the two variations e± 
and e n themselves. 



1. Three phase capillary equilibrium 

The simplest example of a three phase line between de- 
formable surfaces occurs when three capillary interfaces 



meet, for instance at the three-phase-line between three 
mutually immiscible liquids 1, 2, and 3, having mutual 
surface tensions er 12 , a 23 and a 31 . In this case no ad- 
hesion energy is involved (or, alternatively, it may be 
considered as part of the surface tension). The contact 
line variation thus consists of three identical boundary 
variations 

sh c1 = -|d.s{e/ i2a +e/ 23a +e/ 3la }^ 

= J ds [a 12 l 12 + a 23 l 23 + <7 31 1 31 }SX , (32) 



from which we immediately find the boundary condition 



a 31 =0 



(33) 



This expresses nothing but the force balance between the 
three directional line tensions er 12 = a 12 l 12 etc. and is 
known as the Neumann triangle [l[ . The vector equation 
([33]) corresponds to two scalar equations (since there is 
no component along t). These are sufficient to deter- 
mine the three contact angles between the three phases 
(because their sum equals 360°). Notice that this con- 
versely implies that by measuring these angles one can 
only determine the ratios between the three tensions, not 
absolute values. How all this information is conveniently 
extracted is discussed in detail in Ref. [3, Chap. 8]. 



2. Adhesion of two vesicles 

For the case of two adhering vesicles we assume that 
vesicle 1 has bending modulus k\ and tension a%, while 
vesicle 2 has corresponding values k 2 and a 2 . If the two 
bilayers can slide past each other in the region where they 
adhere, their joint bending modulus is given by K12 = 
Ki + «2, because the energies required to bend either 
one just add; the same applies to the tension: o"i2 = 
01 + o~2- We will for simplicity look at the case where the 
spontaneous curvature is zero. The contact line variation 
now contains one adhesion term and three free boundary 
variations. Using the decomposition of SX and Se b as 
given in Eqns. ([30]) and (f3"Tj) . respectively, we find the 
total energy change to be 



SH cA = - d 



-w~ 1 -(k 1+ k 2 ){K 12 ) 2 



l ri {Kl) 2 + 1 -n 2 (Kl) 2 



-Vj 



KiK]_ + K 2 K\ - (Ki + K 2 )K 



12 



K X Ki + k 2 K\ - {ki + K 2 )Kf K ± e 



Kl K{ 



K 2 K\ 



(ki + K 2 )K 12 



(34) 



All corresponding Ku contributions cancel, since K\i is 
again continuous across C; the same happens to the ten- 
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sions. The four terms belonging to the independent vari- 
ations e±, e n , K ± e i , and Vj_e n must vanish individu- 
ally. Since the last two terms have identical prefactors, 
we finally end up with three boundary conditions: 



Kl (Kl) 2 + n 2 (Kl) 2 - ( Kl + k 2 ){K 12 ) 2 
KiK]_ + n 2 Kl - (ki + n 2 )K 12 

Vx (kiK\ + n 2 K\ - (m + n 2 )Kl 2 



2w , (35a) 
, (35b) 

. (35c) 



Notice that three conditions is just what we need in or- 
der to fix three curvatures. Yet, contrary to the case of 
vesicle adhesion to a substrate, Eqn. (|22p. these condi- 
tions also contain one which involves derivatives of the 
curvatures, namely (|35c|) . The origin of this term stems 
from the e„-variation, which is forbidden for the case of 
a rigid substrate. Since the normal variation also multi- 
plies the normal component of the stress tensor, which (as 
Eqn. (|13p informs us) always contains one more deriva- 
tive than the tangential one, this brings about the higher 
derivative condition. 

Regrettably, Eqns. (j35|) do not look particularly 
transparent. It is however possible to symmetrize 
Eqns. (|35al35b[) and thus obtain the two more sugges- 
tive equations 



i + -*3 



2w 
2w 

K 2 



(36a) 
(36b) 



From Eqn. (|35b[) it follows that one of the Kj_ is bigger 
and the other one smaller than K 12 . Hence, when taking 
the square root in Eqns. (|36p . exactly one of the two will 
necessitate a minus sign. 

To conclude this section, let us look at two special 
cases of these boundary conditions which turn out to be 
quite instructive. First, if k 2 — > oo, the second vesi- 
cle approaches the limit of a rigid substrate. In this 
case Eqn. (|36bj) shows that K] 2 = K\ and Eqn. (|3"6li|) 
reduces to the old contact-curvature-condition we had 
just derived for rigid substrates, Eqn. (f2"2")l . And the 
curvature of this effective substrate is determined from 
V \_{K\ — K 12 ) = 0. This latter condition shows that 
the "substrate" -curvature is even differentiable across C 
- or, in other words, the "substrate" shape is a three 
times continuously differentiable function. 

And second, if the two membranes have identical bend- 
ing moduli Ki = k 2 = k, a "symmetrized" contact cur- 
vature condition ensues which reads 



(Ki-Ki 2 )- 



K 2 



K 



12 



W 



(37) 



which tells us that the (squared) curvature jump 2w/k 
demanded by the rigid substrate version ([22| is shared 
between the two membranes, while the final condition 
becomes V±(Kl + K\- 2K] 2 ) = 0. 



V. SUMMARY 

We have shown how the boundary conditions pertain- 
ing to the contact line between a fluid surface adhering 
to a solid substrate or another deformable surface can 
be extracted from a systematic boundary variation in a 
completely parametrization independent way. We would 
like to close with a summary of our main results and some 
remarks: 

• Integrability of the surface energy density Ti. en- 
forces continuity of certain geometric variables 
across the contact line. 

• The highest derivative in Ti thus dictates which 
geometric variables may change discontinuously 
across C in response to adhesion. Hence, for the 
Hclfrich Hamiltonian the tension a does not enter 
the boundary condition even if it enters H; likewise, 
neither tension a nor bending modulus k enter the 
boundary condition if also a gradient-curvature- 
squared term is present in Ti. 

• Higher order derivatives in Ti create boundary 
terms in the variation which pick up surface vari- 
ations that are one order lower. If the curvature 
enters Ti, then a change in slope is noticed, if a 
gradient in curvature enters Ti, then changes in 
curvature are noticed. For this reason the capil- 
lary Hamiltonian is the only one which only picks 
up translations, such that the energy minimiza- 
tion can be reinterpreted as a force balance. In all 
other cases higher derivative deformations (such as 
torques or even more complicated constructs) con- 
tribute to the boundary variation. 

• Generalizations to surfaces hosting additional 
scalar or vector fields (such as composition or tilt 
order) appear straightforward, since these are read- 
ily incorporated into the present framework 



• To be sure, knowing the boundary conditions does 
not mean that one also knows the position of the 
contact line. Rather, the latter has to be deter- 
mined simultaneously with the surface shape. In 
general this task is very difficult, but it is not the 
subject of the present work. 
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